Galaxy Groups and Clusters 


One major difficulty with group/cluster identifications is the con- 
version of velocity to distance. In general, this is done via the 
adoption of a Hubble Constant. However, 


e Peculiar velocities can affect distance measurements by typi- 
cally ^» 4 Mpc 

e Galaxy clusters can have a large velocity dispersion. This cre- 
ates a "finger of God" effect, where the galaxies appear to line 
up along the line of sight. 

e The Hubble flow of large superclusters may just now start to 
reverse. This can cause the galaxies to appear flat against the 
sky. 
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[Huchra et al. 1990, Ap.J. Supp., 72, 433] 


Group/Cluster Definitions 


[Huchra & Geller 1982, Ap.J., 257, 423] 
[Nolthenius & White 1987, MNRAS, 235, 505] 
[Yang et al. 2005, MNRAS, 356, 1293] 


The traditional approach to defining a galaxy group is the per- 
colation, or friends of friends approach. Simply put, any galaxy 
that is within a distance r of another galaxy is a member of that 
galaxy’s group. So if galaxy A is close to galaxy B, and galaxy 
B is close to galaxy C, then A, B, and C are all members of the 
same system. This definition is straightforward, but it depends 
critically on the connection length, r. 


A second method is to define a group as a region where the 
density enhancement of galaxies is a factor of n above the back- 
ground. One can use the friends of friends approach to create 
an initial catalog, and then use the density enhancement criteria 
(over some assumed shape) to define r. 


A completely different method invokes the results of CDM cluster 
simulations, and models of heirarchical cluster formation. One 
first assumes a mass-to-light ratio, and tentatively defines group 
membership (perhaps with a friends of friends approach). One 
then uses the total group luminosity and the mass-to-light ra- 
tio to estimate the group’s mass. Once this is done, the halo 
occupancy predictions of CDM models are used to estimate the 
group's size, and therefore its velocity dispersion. Finally, the de- 
rived velocity dispersion is compared to the observed dispersion, 
and group members are added or subtracted to improve agree- 
ment. The process is then repeated until the group definition 
converges. 





Virgo and the Virgocentric Flow 


[Huchra 1988, ASP Conf. 4, The Extragalactic Distance Scale, 257] 
[Binggeli, Tammann, & Sandage 1987, A.J., 94, 251] 


The dominant system in the Local Supercluster is the Virgo Clus- 
ter. However, it is not a relaxed system it can be divided into 
numerous subsystems. 
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Fic. 4. General outline of the 
8 8 Virgo cluster, based on all clus- 
g ter members listed in the Virgo 
Cluster Catalog. (a) Plot of all 
5 6° i members. (b) Isopleths based 
& : on number counts in cells of 
« : size 075x075, smoothed by 
Z : averaging over 1°5X1°5. The 
a : contours are labeled in units of 
Lu emcee M galaxies per cell, i.e., number of 
Q 16 galaxies per 0.25 square degree. 
(c) Same as (b) but with lumi- 
nosity-weighted isopleths, la- 
14 beled in units of 10'? L per 
deg^The positions of M87, 
M49, and M59, which are asso- 
ciated with the three major 
12 clumps A (M87 cluster in Fig. 
1), B (M49 cluster), and C are 
indicated. For orientation, the 
10 position of the supergalactic 
equator is given in (b). 
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Fic. 23. Velocity distribution for the entire Virgo cluster sam- 
ple broken up into four different morphological groups, and 
into galaxies of all types in the areas of cluster A (around M87) 
and cluster B (around M49), defined by the two contours la- 
beled “6” in Fig. 4(b). The arrows pointing downwards mark 
the corresponding velocity means, which are also listed in Ta- 
ble IV. For comparison, the velocity of M87 is indicated by an 
arrow pointing upwards. 





1500 to 3000 


Figure 8. Surface distribution of (a) all galaxies with veloc- 
ities < 3000 km s~'in the Virgo core, (b) the same for those 
galaxies with velocities less than 500 km s~', (c) galaxies be- 
tween 500 and 1500 km s^!, and (d) galaxies with velocities 
between 1500 and 3000 km s^!. The clumpy nature of the 
cluster core is easily seen. (From Huchra 1985). 
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Fic. 2. Differential luminosity functions for the 
five concentrations appearing in Fig. 1. All 
E + SO + Spiral galaxies within each circle in 
Fig. 1 are counted; dwarf galaxies (as classified 
in the Catalog) are not included. Open histo- 
grams refer to galaxies classified as members in 
Paper II, hatched histograms to possible 
members. 


The Virgocentric Flow 
[Aaronson et al. 1982, Ap.J., 258, 64] 


The infall of galaxies towards Virgo dominates galactic motions 
in the Local Supercluster. As illustrated by this apple diagram, 
in some directions, the radial velocities of galaxies are tripled- 


valued. 
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Fic. 3.—Normalized redshift residuals for solution 1, roughly log Uprea — log vos, with positive residuals pointing away from the position of 
the Local Group at z — 0. Galaxies are plotted on a meridional plane, with positions inferred from H 1 velocity widths. Contours show 
predicted redshifts, uncorrected for the peculiar motion of the Local Group, spaced at intervals 1/10 the Virgo redshift (after correction for 


our peculiar z velocity). 





Large Scale Bulk Flows 
[Tonry et al. 2000, Ap.J., 530, 625] 


The Local Group is moving towards Virgo (bł! = 74?,1!! = 283?) 
at a velocity of ~ 270 km s^ !; this is called the “Virgocentric 
Infall.” However, when one measures the anisotropy in the mi- 
crowave dipole radiation, the absolute motion of the Milky Way 
is ~ 620 km s^ in the direction of b/! = 27°, £ = 268°. This 
is presumably due to the potential of the Hydro-Centaurus Su- 
percluster (once known as the Great Attractor). Unfortunately, 
the center of this supercluster lies in the Galactic plane, so it is 
difficult to study. 











80 ~r [T3 


(Mpc) 


"SQ 
































—80 —60 


"SGX" (Mpc) 


Fic. 21.—Same as Fig. 20, except that our survey galaxies are overplotted. Note that the Centaurus galaxies that appear to be going through the Great 
Attractor actually pass above it by about 15°, and hence, lie in the stall zone. 
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The Two-Point Correlation Function 


The most common way of describing the overall clustering of 
galaxies is through the two-point correlation function. This can 
be performed either on the projected distribution of galaxies on 
the sky (the 2-D case) or the 3-D galaxy distribution. 


The premise of the two-point correlation function is simple. If 
galaxies are Poissonly distributed, the number of galaxies you 
can expect to find in a volume dV; is 





where po is the average density of galaxies in the entire volume. 
(If dV, is infinitessimal, then this number is a probability, since 
the chance of finding more than one galaxy in the volume is zero.) 
Similarly, the probability of finding a galaxy in volume dV; and 
finding another galaxy in volume dV» is 


dNpair = pgdVidVo (23.02) 


We can now define the two-point correlation function, €, as the 
probability of actually observing such a pair of galaxies, com- 
pared to the Poisson probability, t.e., 


dNpair = poll + €]dVidV2 (23.03) 


where € is a function of the separation between the two volume 
elements, r. In other words, if €(r) > 0, the positions of the 
galaxies are correlated; while if £(r) < 0, the positions are anti- 
correlated. 


Let’s look at £(r) another way. Let ó(x) be the relative overden- 
sity of region x, i.e., 


5(x) = LL (23.04) 


The probabilty of finding a pair of galaxies in volumes dV; and 
dV5 centered at positions xı and zə is 


dN pair = p(xı)dVı ] p(z2)dV» 
= po {1 + ó(z1)) - po 11 + 9(22)) dVı dV 
= po {1 + 0(21)) {1 + 9(23)) dVi dV2 


= po {1 + 6(x1) + 6(x2) + 0(21)6(x2)) dV, dV5 
(23.05) 





Averaged over the entire field, ó(z) must equal zero, so 
iN pair = pé {1 + 6(x1)6(a2)} dV, dV5 (23.06) 


If r is the separation between position x; and position x2, then 
a comparison of (23.03) and (23.06) yields 


&(r) = (6(z)0(z  r)) = (87) (23.07) 


In other words, £(r) is a measure of the overdensity as a function 
of radius. 
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Computing the Two-Point Correlation Function 





In order to compute the 
two-point correlation func- 
tion, one has to count the 
number of sources between r 
and r + ôr from every object 
in the catalog, and compare 
that number to the number 
you would get from a random 
distribution. 














This latter step is tricky. It is usually impossible to estimate the 
signal of a random distribution analytically: bright stars, changes 
in the detection threshold (perhaps from extinction), and other 
factors will effect the ^window function" of the observations. In 
addition, most surveys have edge effects, and objects close to the 
boundary have fewer neighbors than those near the center. While 
one can try to compute the correlation function using subsets of 
data far from the survey boundaries, or do a correction for in- 
complete shells (Rivolo 1986), most investigators simply perform 
Monte Carlo experiments. One randomly adds points to the sur- 
vey area using the exact window function as the original data, 
and computes 





N,a DD(r) 
O= N RRO 








SES (23.08) 


where N is the number of data points in the catalog, Nya is the 
number of random points, DD(r) is the number of pairs with 
separation r in the data catalog, and RR(r) is the number of 
pairs with separation r in the random data set. (To minimize 
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the errors of this procedure, the number of random data points 
must be many, many times that of the data.) 


Unfortunately, this type of procedure still does not remove all the 
edge effects from the calculation. To remove this bias, a series of 
improved estimators have been proposed, including one by Davis 
& Peebles (1983) 








e(r) = EL k Ed (23.09) 
Hamilton (1993) 
o eee MM AM a (23.10) 


and (currently the most popular) Landy & Szalay (1993) 








RR(r) N RR(r) 








2d] DD(r) , Nra DR(r) (23.11) 


pE 


where DR(r) is the distribution of pair separations in a data 
versus random point catalog. All these estimators have some 
bias to them, as they assume a mean density that is defined from 
the sample data, rather the true, universal mean. 
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Angular versus Spatial Correlation Functions 


Until recently, large extragalactic datasets with full 3-D informa- 
tion (i.e., with redshifts) were rare. Thus, most early analyses 
involved the angular correlation function, w. One can infer the 
three dimensional correlation function from the 2-D function with 
a little mathematics. Suppose the true 3-D radial distribution of 
a galaxy cluster with radius Rpg, can be written in the form 


n(r) = no [1 + £(r)| dV (23.12) 


> 





Then, projected on the sky, this distribution will be 


N(R) =2 J us agar =2 J IA d 
R (r2 — R2) /2 R (r2 — gy? 
(23.13) 
where R is the projected distance from the center of the cluster. 
Now suppose the correlation function can be described as a power 
law, i.e., £(r) x r ?. Then 


zuo 


N(R) =2 J an [1 +ar]r (r? -Rdr — (2314) 


The first term of this expression integrates to a constant offset. 

The second term is complicated, but if Rmax —> oo, then 
N(R) x R^ (23:15) 

Thus, a spatial correlation function which is a power law will 


produce an angular correlation function that is also a power law, 
with an exponent that differs by one. 


13 


Results from the Two-Point Correlation Function 


The results of various galaxy surveys show that the two-point 
correlation function for galaxies is approximately a power law in 
the range 0.1 Mpc « r « 16 Mpc, with 


Sirve n fng (23.16) 


and y ~ 1.75. The value of ro is called the correlation length. For 
r < ro, (r) is positive, hence galaxies are clustered on smaller 
scales. For r > ro, €(r) < 0, implying a deficit of galaxies in 
this range. Observations indicate that rg ~ 6h~' Mpc. Inter- 
estingly, the correlation function for galaxies, galaxy groups, and 
rich clusters all have approximately the same power-law index. 


The only difference is their correlation lengths. 
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Figure 4: Galaxy 2-point correlation function at the present epoch. Red symbols (with vanishingly 
small Poisson error-bars) show measurements for model galaxies brighter than My — —23. Data for the 
large spectroscopic redshift survey 2dFGRS?? are shown as blue diamonds. The SDSS?^ and APM?! 
surveys give similar results. Both, for the observational data and for the simulated galaxies, the corre- 
lation function is very close to a power-law for r < 20/1-! Mpc. By contrast the correlation function for 
the dark matter (dashed line) deviates strongly from a power-law. 


[Springel et al. 2005, Nature, 435, 629] 
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The Correlation Correlation in Fourier Space 


Rather than work in terms of £, one frequently wants to know 
how much clustering there is at a given linear scale length. This is 
given by the Fourier transform of the correlation function, P(k), 
where k is the wavenumber of the fluctuation (i.e., 27/A). The 
two are related via 








Y qoos 
P(k) = > / g(r) m 4rr?dr (23.17) 
and T TI 
= sin Kr 2 
Er) - Z3 J P(k) =k dk (23.18) 


In this form, we can express the correlation function and density 
fluctuation in terms of mass. First, we note that if the initial 
spectrum of fluctuations had no preferred scale, then the simpli- 
est form for the power spectrum would be a power law, 


PR ook” (23.19) 
If we put this in (23.18), then 
£(r) x n EET EE (23.20) 
0 kr 


We can quickly estimate this integral by noting that when kr > 
1, the integrand rapidly goes to zero. Moreover, when kr < 1, 
sin kr/kr ~ 1. So to first order 


1/r 


Lr 
E(r) ox i kdk ak ocr ta) (23.21) 
0 





0 


Since M «x pr?, 
£(M) x M7("+3)/3 (23.22) 
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Moreover, since € has already been associated with the root mean 
square of a density fluctuation via (23.07), the mass spectrum of 
density perturbations should be 

(83,) ^ = AM- (1/6 (23.23) 


where A is just some constant of proportionality. 
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Growth of Structure 


At the time of the Big Bang, the universe was extremely smooth, 
but today it is extremely clumpy. How did we get to this state? 


First, let's define two classes of perturbations. Isothermal per- 
turbations are fluctuations that affect the matter of the universe, 
but not the surrounding radiation field. These are the types of 
fluctuations that occur today, since matter and radiation are de- 
coupled. But in the early universe (in the era before recombi- 
nation), this type of fluctuation may not have existed. 





The second class of perturbations are called adiabatic perturba- 
tions. These effect both light and matter together. Such per- 
turbations do not occur today, but they might have occurred at 
early times before decoupling. 
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Isothermal Perturbations 


Fluctuations grow due to gravity; the gravitational force per unit 
mass is 


e = GPR (23.24) 


Now consider an isothermal perturbation, where gravity is at- 
tempting to move electrons (and protons) through a radition 
field. The velocity of the electrons is v; the cross section of the 
electrons to photon absorption is or (the Thomson cross sec- 
tion). Because the electrons are moving through a medium filled 
with photons, they will experience a drag force. Since the pro- 
tons must follow the electrons, they will also (indirectly) feel this 
force. Per unit mass, the drag force is 


I xs 





aT^ (orv) 


23:95 
mos (23.25) 


Peas a 


(It’s not difficult to see where this equation comes from: aT“ 
is the energy density (ergs/cm?) of the microwave background, 
and ov is the cylindrical volume swept out by the electron each 
second. So the numerator is the amount of energy encountered 
by the electron each second as it plows through the photons.) 


--------- Vt---------- 
O 
Now, since the universe is expanding and photons are being red- 


shifted, 
T=T%(1+z) p=po(l+z)? (23.26) 
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If we compare Farag to Foray, and explicity put in the these de- 
pendencies, then 


Farag m 3oTvaT^ u 3orvaT$(1 ee p (23 2n) 
Foras i 4r Gm gcRp B An Gmg cRpo(1 + 2)» : 





Now consider R, the physical size of the perturbation being con- 
sidered. In an Einstein de-Sitter universe, the time it takes an 
electron to move a distance R is approximately vt. If perturba- 
tions begin at the time of Big Bang, then by redshift z, the size 


1S 
2 1 
= wd 3? 23.2 
R 3H, + z) (23.28) 


If we use this to substitute for R, then 
Frag 9a aT Ho 


Een 9m 23.29 
Po STEM mcpe fo Up) 





Numerically, this works out to be 


Fara z 
N bes) (23.30) 


grav 


Decoupling occurred at z ~ 1500. Before this time, Farag > 
Foray, 80 gravity could not overcome the drag force of the thermal 
background. So isothermal perturbations did not grow during 
this era. 
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Adiabatic Perturbations 


If the radiation field is perturbed along with the matter, the ex- 
cess energy contained in a region will diffuse outward and damp 
out the fluctuation. However, this damping does not occur in- 
stantaneously. 


Let A = 1/(neor) be the mean free path of a photon. In terms 
of A, the time between scatterings for photons is A/c, and, from 
statistics, the number of scatterings necessary for a photon to 
random walk a distance R isa (R/A)?. For a fluctuation not to 
be damped out, the diffusion time must be longer than the age 
of the universe, i.e., 


RN TAN 21 E. WU -3/2 

ch qe suc — 2 

(x) (322042 doa FA 
(23.31) 


Now if we substitute mass for radius and again note that the 
protons must follow the electrons, so that p = nemy, then the 
condition for perturbation growth is 


2cmy 


3H por 


3/2 
: 5 (1+ d (23.32) 


4 4 
M = =1pR? > sm 
or, if we use (23.26) to substitute for density, 


4r [ 2emg | ^? pp 15 
— Legare q983 
M > 3 tae} 2:06 (1-4 z) (23.33) 








where pe is the critical density of the universe. At decoupling, 
this works out to M > 10!°Mo. This is somewhat of an un- 
derestimate, since A is changing rapidly in the early universe. 
But it does demonstrate that adiabatic perturbations can only 
propagate if they are very large; small scale perturbations will 
be damped out. 
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The Jeans Mass 


Another question we can ask concerns the minimum mass for 
gravitational collapse. For collapse to occur, the gravitational 
potential energy must overcome the thermal motion of the gas. 
If v, is the sound speed, 


G 
ie 


7 v? (23.34) 


N| = 


If we substitute density for radius, then 


iron 1/3 : 
2 ——— 23. 
om (4) > v? (23.35) 
SO : 
E ed 
Mz = (=) G p. (23.36) 


where My is called the Jeans mass. After decoupling, matter 
acted as an ideal gas, so 


v= Can - CM (23.37) 


p MH 


If we plug in the numberes for the universe shortly after decou- 
pling, i.e., z ~ 1000, then v, ~ 4 km s^! , and the Jeans mass is 
M; ^ 10°Mo. This is similar to the mass of a globular cluster. 


In the early universe, however, the ideal gas law did not apply. 
At z ~ 40,000, radiation pressure dominated, and the energy 
density was much greater than the matter density. So 


T ps 
plc. ang p= aa (23.38) 


c? 
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During this time, the sound speed was v, ~ c/V/3, and the density 


Was fi 
p= Te 3-2)! (23.39) 


The Jeans mass was therefore Mi; ~ 10!9.M(5. Again, this 


demonstrates that large fluctuations could grow in the early uni- 
verse. 
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Growth of Perturbations 


Can we connect the structure we see today to perturbations in 
the microwave background? The largest structures in today’s 
universe are superclusters; these have overdensities of óp/p ~ 2. 
To see what that means, consider a region of slightly enhanced 
density, dp, within an Einstein de-Sitter universe. Recall from 
our discussion of the basic equations of cosmology (2.02), the 
expansion of any region of space is given by 


AUD 
R 8 kc? Ac? 
NERO: reete oii 29.4 
Bü 37 G(p9 + 9p) R2 + 3 (23.40) 
which, for p > pe, has the solution, 
R-—a(1—cos0) t= aC) — sin 0) (23.41) 
C 


Let's just consider the R part of the solution. In the early uni- 
verse, 0 was small, so we can Taylor expand the trignometric 
terms as 


1 1 
cos = 1 — 20 4 ri +... (23.42) 
So (23.41) becomes 
Real P ee (23.43) 
uid 51 ve 


If we assume the perturbation is small and neglect the higher 
order terms, then 


= a eT 23.44 
R+0R a (50 5 s) (23.44) 
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SO 


Ro«60? and R x 0* => oF x x R (23.45) 
Now let's relate this to the density. 
M 3M 3p 
m óp „ÖR 
p 
— = 3— 23.47 
r R (23.47) 
õp ôR 
7 x OX Ræ (1+2)! (23.48) 


In other words, the amplitude of a small density perturbation 
will grow linearly with the size of the universe due to the Hubble 
expansion. This puts a limit on the amplitude of the initial 
density perturbations: for example, if a supercluster now has 
0p/ p ~ 2, at decoupling, the initial density fluctuation must have 
been at least 1000 times smaller. 


Note, however, that very early on, when the universe was dom- 
inated by radiation, so p c R^ (since the dominate source of 
matter ter/energy was being redshifted). So at this time 


a R (23.49) 
p 





One can relate these to time, t in the matter dominated era by 
(1.29) 

Ó 

IP x Roc t? (23.50) 
p 
At earlier times, when radiation pressure was important, R went 
as t!/? (that’s an exercise to the student), and 


» x R^ xt (23.51) 
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Press-Schechter Mass Function 


[Press & Schechter 1974, Ap. J., 187, 425] 
[Longair 2008, Galaxy Formation] 


According to CDM scenarios, galaxies and clusters were built via 
hierarchical clustering. A standard analytical treatment of this 
process was presented by Press & Schechter in 1974. Although 
their analysis assumed a critical density universe with no cosmo- 
logical constant, no dark matter, and neglected a host of very 
important effects, its result is very close to that predicted by the 
largest, most modern N-body simulations. 


Assume that the primordial density perturbations were randomly 
distributed in a Gaussian manner. There were just as many 
positive fluctuations as negative fluctuations, and the probability 
of a given fluctuation with total mass M being over (or under) 
dense by an amplitude ôm = (p — po)/po was 


1 ô? 
plm) = Y exp | | (23.52) 


Obviously, since the fluctuations are Gaussian, the mean-square 
fluctuation 


(654) — o? (Mt) (23.53) 


(This is the same description as for the correlation function, ex- 
cept now we are studying Ó as a function of mass, rather than 
distance.) It is also obvious that before collapse, the volume en- 
compassing this perturbation was related to the mean density of 
the background by M = po V. Now assume that once a pertur- 
bation develops an amplitude greater than a critical amplitude 
dc, it evolves rapidly into a bound structure of mass M. (This 
conveniently neglects all the messy physics of the situation.) For 
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a fluctuation with a given mass, the fraction of objects that will 
become bound at a particular epoch is 


1 ee 6° 1 
PM) = roa h, °° [arm] = a een 


where t = ó,/V/20 (.M) and erf is the error function. Now recall 
that the mass spectrum of density perturbations should be 





(844) ^ = AM- 9e (23.23) 
so, if we combine this with (23.53), we get 
c(M) = (84) ^ = AV 4-36 (23.55) 
So 
dig oe M(3tn)/6 — m (23.56) 
V2e(M) | V2A17 : 


where we have replaced ôe with a reference mass .M* 


3/(3+7n) 
M* = (=) (23.57) 


Ue 
The space density of these perturbations is given by 
1 d F(.M) po dF(M) 
= PDT ERE 
Mes 


(The minus sign comes from our convention — the density of 
bound structures is inversely dependent on their mass.) Solv- 
ing for this using (23.54) yields 








N(M)dM = dM (23.58) 


(3n) (3+7) 
3 


NU = gr (145) ae (ae) E. “(aa 








(23.59) 
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All we need now is an estimate of how M* changes with time (or 
redshift). For an Einstein-de Sitter universe, perturbations grow 
as the size of the universe, R, so 


e x Ro t? — (62) =07(M) = AM- C08 & 1/5 


(23.60) 
Therefore 
AA Bea S HL p oam (23.61) 
OT 
" 4/(34-n) 
M* = Mo (=) (23.62) 





where Mo is the value of M* at the present epoch, to. Thus, for 
a given n, we have a prescription for the number of bound objects 
of a given mass versus redshift. For reference, the canonical value 
of n is that proposed by Harrison (1970) and Zel’dovich (1972), 
n= 1. 





There are many, many problem with the Press-Schechter formu- 
lation (too many to describe here). It neglects peak superpo- 
sitions (where one falls on top of another), assumes all objects 
are spherical, and ignores the fact that once structures begins to 
form, non-linear dynamics will begin to dominate. Nevertheless, 
it does a remarkably good job a number of observations. For 
example, note the form of the function: at the high-mass end, 
there is an exponential cutoff (with the power of M/M* close to 
one), at low masses, the distribution predicts a power law. This 
is exactly the form of the Schechter luminosity function. 
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Figure 2: Differential halo number density as a function of mass and epoch. The function n(M, z) gives 
the comoving number density of halos less massive than M. We plot it as the halo multiplicity function 
M?p-! dn /dM, where p is the mean density of the universe. Groups of particles were found using 
a friends-of-friends algorithm with linking length equal to 0.2 of the mean particle separation. The 
fraction of mass bound to halos of more than 20 particles (vertical dotted line) grows from 6.42 x 1074 
at z = 10.07 to 0.496 at z = 0. Solid lines are predictions from an analytic fitting function proposed in 


previous work! !, while the dashed lines give the Press-Schechter model!* at z = 10.07 and z = 0. 


Millenium simulation versus Press-Schechter formulation 


[Springel et al. 2005, Nature, 435, 629] 
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Fig. 16.5 The evolution of the comoving number density of dark matter haloes with 
masses greater than M as a function of redshift for a standard Cold Dark Matter 
model with No = 1. The curves have been derived using the Press-Schechter form 
of evolution of the mass spectrum which is a good fit to the results of N-body 
simulations. The dotted line labelled ó* shows the present number density of L* 
galaxies. (after Efstathiou 1995). 
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